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. Vacuum of Quantum Chromodynamics in very strong (hadron-scale) mag- 

VO \ netic fields exhibits many interesting nonperturbative effects. Some of these 

effects can be studied with the help of lattice simulations in quenched QCD. 
We review our recent results demonstrating that very strong external magnetic 
fields lead to (1) the enhancement of the chiral symmetry breaking [the quark 
condensate rises with the increase of the external magnetic field]; (2) the chiral 
magnetization of the QCD vacuum [spins of the quarks turn parallel to the ex- 
i ternal field]; (3) the chiral magnetic effect [a CP-odd generation of the electric 

current of quarks directed along the magnetic field]; (4) a CP-odd generation 
■ of the anomalous quark electric dipole moment along the axis of magnetic field. 

^ , The first three effects were already predicted theoretically, and subsequently 

I observed numerically in our simulations, while the fourth effect is a new result. 

00 ' 

0^ ■ 1 Introduction 

o 

O ■ Motivation 

Noncentral heavy-ion collisions may create very strong magnetic field due to 
^ ' relative motion of electrically charged ions and the products of the collision. 

This fact ignited the interest of the scientific community towards an inves- 
tigation of the properties of the strongly interacting matter exposed to the 
magnetic field. The QCD effects should be visible because the strength of the 
created magnetic field may be of the order of the hadron-scale, or higher. For 
example, at first moments (r ~ Ifm/c) of Au-Au collision at the Relativis- 
tic Heavy Ion Collider (RHIC) the strength of the magnetic field may reach 
eB ~ (10 — 100 MeV)'^ The strong magnetic fields will also be created in 

the future experiments at the Facility for Antiproton and Ion Research (FAIR) 
at GSI, at the ALICE experiment at LHC, and at the Nuclotron-based Ion 
CoUider fAcility (NICK) at Dubna. 



*Talk given by M.N. Chernodub at the 10th Workshop on Non-Perturbative QCD, June 8-12, 
2009, Paris, France. 
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Vacuum effects due to strong magnetic field in QCD 

There are various vacuum effects wfiicli appear due to the strong external mag- 
netic fields (here we do not discuss interesting phenomena that may emerge in 
a dense quark matter). 

Enhancement and shift of the chiral phase transition. Theoretically, 
the very strong magnetic fields may significantly modify the QCD phase dia- 
gram 131 [I]- The external magnetic field increases the transition temperature 
and enhances the strength of the (phase) transition from the chirally broken 
(low temperature) phase to the chirally restored (high temperature) phase. 
Moreover, in the external magnetic field the transition - which is a smooth 
crossover in the absence of the fields - becomes a first order transition plHj. 

Enhancement of the chiral symmetry breaking. The magnetic field 
stabilizes the chirally broken (low-temperature) phase of QCD by enlarging the 
value of the chiral condensate and enhancing the chiral symmetry breaking. 
This result was obtained in various approaches using the chiral perturbation 
theory [SlEllTj, the Nambu-Jona-Lasinio model [8j, the linear sigma model [9] 
and the AdS/QCD dual description [10]. In lattice simulations this effect was 
observed in our recent work [11] . In Section [2] we discuss main results of [TT] . 

Chiral magnetization. Another effect of the strong enough magnetic 
field is a (para)magnetic response of the QCD vacuum: the external fields po- 
larize the spins (or, equivalently, the magnetic moments) of the quarks and 
antiquarks, leading to the chiral magnetization of the vacuum. At relatively 
low magnetic fields the chiral magnetization is proportional to the strength 
of the magnetic field. The coefficient of the proportionality, called the chi- 
ral magnetic susceptibility, was first discussed in Ref. [I2] in order to ana- 
lyze phenomenologically interesting nucleon magnetic moments. The value of 
the magnetic susceptibility - which was estimated using various analytical ap- 
proaches [131 [HI [151 [T71 [18] ~ can be measured in experiments on lepton pair 
photoproduction [19], in radiative heavy meson decays [20], etc. At higher mag- 
netic fields the behavior of the magnetization deviates from a linear function as 
it gets affected by logarithmic corrections at moderately strong fields [21]. At 
asymptotically strong fields the magnetization reaches an upper bound. The 
first lattice investigation of the chiral magnetization was performed in Ref. [22] . 
The results of [22] are briefly described in Section [3l 

Chiral Magnetic Effect. The magnetic fields may lead to quite unusual 
effects due to the nontrivial topological structure of the QCD vacuum [H |23] . 
We discuss a particular realization of the Chiral Magnetic Effect (CME), which 
gives rise to a generation of an electric current along the direction of the mag- 
netic field in a nontrivial topological backgrounds of gluons [Tj. The potentially 
observed feature of the CME is a non-statistical asymmetry in the number of 
positively and negatively charged particles emitted on different sides of the re- 
action plane in the heavy-ion collisions [23]. There are preliminary indications 
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that this CP-odd effect has been indeed observed by the STAR Collaboration 
in experiments at RHIC |25j . In our lattice studies we found the existence of 
the CME both at vacuum and thermal ground states of gluon fields, and at 
specially prepared configurations with nontrivial topological charge |26j . The 
main results of Ref. [26] are given in Section [H 

Quark electric dipole moment. Another manifestation of a nontrivial 
topological structure of the QCD vacuum is an appearance of an electric dipole 
moment of a quark along the direction of the external magnetic field. This CP- 
odd effect - which is a spin analogue of the Chiral Magnetic Effect - is observed 
in lattice simulations in Ref. [27J. The lattice evidence of the generation of the 
quark's electric dipole moment is discussed in Section [5l 

Physical setup and technical details of simulations 

We utilize lattice QCD with the simplest SU (2) gauge group because all studied 
effects originate in the chiral sector of QCD where the number of colors is not 
crucial. In our simulations only valence quarks interact with the electromag- 
netic field, and the effects of the virtual quarks on gluons are neglected. Indeed, 
the inclusion of dynamical (sea) quarks makes the simulations computationally 
difficult, while the essential features of all mentioned effects are visible in the 
quenched limit as well (an extended discussion is given in |1H [22l [26l [28] ) . 

The chiral effects are best studied with massless fermions. In order to imple- 
ment chirally symmetric massless fermions on the lattice, we use Neuberger's 
overlap Dirac operator [29j. We reduce the ultraviolet lattice artifacts using the 
tadpole-improved Symanzik action for the gluons. We introduce the magnetic 
field = B6^3 into the Dirac operator by substituting su (2)-valued vector 
potential with u (2)-valued field. In infinite volume such substitution is 
A^fi A^fi + Cfj,5^^ , where Cfj_ = B {x25fj_i — xi5^2) /2 is the external U (1) gauge 
field. In the finite volume with periodic boundary conditions an additional x- 
dependent boundary twist for the fermion fields should be introduced [30], and 
the uniform magnetic field is forced to take quantized values [31j, qB = 2it k/L'^, 
where k £ 1^ and q = 1/3 |e| is the absolute value of the electric charge of the 
d-quark. Note that in all figures below we indicate the magnetic field strength 
in units of qB = eB/3 and not in eB. In our simulations the lowest magnetic 
field takes values V eB ~ (400 — 600 MeV)^. Such fields are stronger than those 
expected at RHIC, while they are of same order as the magnetic fields that will 
presumably be created at the LHC collisions [2j . 

In order to check the finite- volume effects we use two spatial lattice volumes, 
14^ and 16^. We study the system at three different temperatures, T/T^ = 
0, 0.82, 1.12 (the critical temperature in SU{2) gauge theory is Tc ~ 310 MeV). 
In order to make sure that the influence of the finite ultraviolet cutoff is under 
proper control, we utilize three lattice spacings, a = 0.089 fm, 0.103 fm, 0.128 fm. 
Our spatial volumes range from = (1.4 fm)^ to (2fm)^. 
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2 Enhancement of chiral condensate 



The chiral condensate 

S=-(0|gg|0), (1) 

is commonly used order parameter for the chiral symmetry breaking in the 
theory with massless quarks. The condensate is zero if the chiral symmetry is 
unbroken. If the strength of magnetic field is larger than the pion mass (which 
is zero for massless quarks) but is still much smaller than the hadronic scale, 
one can use the chiral perturbation theory to calculate the field dependence of 
the chiral condensate. According to the original work [5], in a leading order in 
B, the chiral condensate S (B) rises linearly with the field strength: 

where F.,^ ~ 130 MeV is the pion decay constant. The corrections due to non- 
vanishing pion mass may also be calculated 0. 

The condensate ([T]) can be calculated using the Banks-Casher formula |32j . 

S = lim hm ^ (3) 

where V is the total four-volume of Euclidean space-time. The density p (A) of 
eigenvalues An of the Dirac operator V = {d^ — iA^) is defined by 

2?^n = A„Vn, p(A) = (^5(A-A„)), (4) 

n 

where ipn are the eigenmodes of the Dirac operator. 

According to Eq. ([3]) the enhancement of the chiral condensate in the mag- 
netic field means that the Dirac eigenvalues in the vicinity of A = should be- 
come denser as the magnetic fields increases. In order to check this fact we plot 
in Figure [1] (left) the dependence of lowest Dirac eigenvalues on the strength 
of the magnetic field. One can clearly see that the increase of the strength 
of the magnetic field leads to the consolidation of the near-zero modes. This 
configuration has topological charge equal to one, and therefore the configura- 
tion contains one exact zero mode, A = 0, in agreement with the Atyah-Singer 
theorem. The presence of the zero mode is independent of the strength of the 
external Abelian magnetic field because the number of zero eigenmodes is equal 
to the topological charge of the gauge field configuration, while the topological 
charge is not affected by the Abelian magnetic field. 

In Figure [T] (right) we show the chiral condensate as the function of the 
magnetic field qB at zero temperature and at T = 0.82 Tc. Following ([2]) we fit 
the chiral condensate by the linear function, 

S(i?) = So + , (5) 
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Figure 1: (left) The spectrum of lowest twelve eigenvalues A (including the zero mode 
A = 0) of the Dirac operator in a background of a typical gluon configuration vs the 
magnetic field qB. (right) The chiral condensate vs qB at two different temperatu- 
res T. The solid lines are the linear fits by the function ([5]). Figures are from 



where Eq and are the fitting parameters. 

The best fits are shown in Figure [1] (right) by the sohd lines. The best fit 
parameters at zero and nonzero temperatures are, respectively [llj: 

r = 0: = [(320 ± 5) MeV]^ = (1.53 ± 0.11) GeV,(6) 

T = 0.82Tc: = [(291 ± 1) MeV]^ Af = (1.74 ± 0.03) GeV. (7) 

We see that the increase of temperature leads to a decrease the chiral condensate 
(as expected), and to a decrease the slope of the S-dependence. 

The zero-field zero-temperature best fit value ([6]) of the chiral condensate, 
Sq*, agrees very well with other numerical estimations in quenched SU{2) gauge 
theory [33]. Surprisingly, the numerical value of the slope parameter Ag is quite 
close to the T = result of the chiral perturbation theory ([2]), Ref. [5]: 

A*^''"'^^ = Ai^F^/^hPi = 1.97 GeV . (8) 

The similarity between the theoretical prediction ^ and first principle T = 
result ([6]) should be taken with care. First, we are studying the quenched 
theory, in which the virtual charged pions are absent, and what we observe is 
the effect of the gauge fields. Second, the strength of our weakest magnetic 
fields is still greater than the scale imposed by the zero-field chiral condensate, 
Eq . In this regime the prediction ([2]) of Ref. |5] should not work, in general, 
as the linear behavior is expected to be realized for much weaker fields. Thus, 
the observed linear enhancement of the chiral condensate in the absence of the 
pion loops is an intriguing unexpected feature of the non-Abelian gauge theory 
exposed to very intense external magnetic fields. 
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3 Magnetization of QCD vacuum 



Quarks are spin-1/2 particles carrying magnetic moments. The magnetic mo- 
ments are polarized by the external magnetic field. A natural quantitative 
measure of the polarization is given by the expectation valu€0 [E] 

(^S,^^) = x{F) (^^) qF^p , (9) 

where Sq,/3 = ^[707/3—7/370] is the relativistic spin operator, = 9^ajy — 9jya^ 
is the strength tensor of the electromagnetic field a^. 

The right hand side of Eq. Q is proportional to the electromagnetic field 
strength tensor due to the Lorenz covariance. The proportionality to the chiral 
condensate {^^) (evaluated at the nonzero external electromagnetic field F) 
allows us to disentangle nonlinear effects of the enhancement of the chiral con- 
densate in the external magnetic field (discussed in the previous Section) from 
the effects of the quark's spin polarization. The strength of the vacuum polar- 
ization is characterized by the chiral magnetic susceptibility xiF)- 

The (chiral) magnetization of the QCD vacuum in the external magnetic 
field B = F12 = —F21 can be described by the dimensionless quantity 

li{qB)=x{qB)qB ^ (^Si2^) = ;u(gS)(^^) . (10) 

In [22j we derived a magnetization analogue of the Banks-Casher formula ([3]): 

(^S„/3^) = / d*x^{{x) ^o^pM^)) ■ (11) 

Moreover, we have proven theoretically in ^22j and also checked numerically, 
that the following factorization rule holds: 

(¥E,/3*) = d\ilj{{x)T.o^piljx{x)). (12) 

The comparison of this formula with Eqs. ([9]) and (jlOp gives 

li{qB) = x{qF)qF^p = lim( j d^x Vl(x; S) S12 Va(x; i?)) , (13) 

where ipx{x]F) is the eigenmode of the Dirac operator in the external (mag- 
netic) background field B = F12. 

We established the factorization property (|12p for all studied values of the 
magnetic fields. The comparison of the left and right sides of Eq. ([12]) evaluated 
at the same set of configurations is shown in Figured (left). While both esti- 
mations agree with each other within error bars, the factorized definition (|13p 
is much more accurate compared to the nonfactorized one (jlip . 

Flavor and spinor indices are omitted in Below we consider one quark flavor for simplicity. 
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Figure 2: The magnetization /i vs the magnetic field qB. (left) Check of the factor- 
ization fll2p : the empty squares and the full circles show the nonfactorized (1111) and 
factorized f|T3|) definitions, respectively, (right) The magnetization /i at two temper- 
atures. The lines correspond to the best fits f|T4l) . Figures are from Ref. [22j . 



We show in Figure [2] (right) the magnetization for two different tempera- 
tures, r = and T = 0.82 Tc. The behavior of the magnetization is consistent 
with general expectations: at low magnetic fields the magnetization is a linear 
function, which indicates the existence of a nonzero susceptibility at vanishingly 
small external magnetic field. At high magnetic fields the quarks are fully po- 
larized and the magnetization approaches the saturation regime, fJ-{qB) — > —1. 

We fit our data by the function (other fitting functions are discussed in |22] ) 

f^(B) = ^J^.,,t.n'-^. (14) 

which has two fitting parameters: the zero-field susceptibility xo and the infinite- 
field saturation parameter /loo- The fits - shown by the lines in Figure [2] (right) 
- give the following values for the chiral susceptibility: 

_ r -1.547(6) GeV-2 Auv~2GeV T = , . 

-1.53(3) GeV-2 Auv~l-5GeV T = 0.82Tc ^ ^ 

where T is the temperature, and Auv is the momentum scale, at which the sus- 
ceptibility is evaluated (the momentum dependence is discussed, e.g., in |14j). 

An experimentally relevant and phenomenologically interesting quantity is 
the product x(^^) at vanishing field |T^. Our data give the following result: 

X {^^) = -46(3) MeV [quenched Hmit] . (16) 

which is surprisingly close to the estimation based on the QCD sum rules tech- 
niques, xC*"^) ~ -50 MeV pL3j. 
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4 Chiral Magnetic Effect 



In brief, the physical mechanism behind the Chiral Magnetic Effect (CME) is 
as follows [Ij. A strong magnetic field forces the magnetic moment of quarks 
to turn parallel to the direction of the field. If the quarks are light (mass- 
less) then the left-handed quarks will move, say, towards the direction of the 
magnetic field while the right-handed quarks will move backwards. If there is 
an imbalance between left-handed and right-handed quarks, then a net elec- 
tric current (x) = tp{x)'y^ip{x) appears along the magnetic field [U [231 El]- 
This longitudinal electric current may lead also to a spatial separation of the 
electric charges. The longitudinal current and the spatial charge separation 
are the consequences of the CME. The chiral imbalance may be created by 
topologically nontrivial configurations of gluon fields. Generally, the CME is a 
reflection of the CT^-odd structure of the vacuum. The existence of the CME 
was demonstrated in lattice simulations in Ref. [26]. 

In Figure [3] we plot a spatial excess of the charge density originating due to 
the applied magnetic field, jo {x; B) = {Jq (x)) b — {jo (x)) b=o- The subtraction 
of the zero-field density removes all ultraviolet contributions providing us with 
a nonperturbative quantity. The magnetic field applied to a typical gluon con- 
figuration leads to a spatial separation of the electric charges in an agreement 
with the CME. The effect increases as the magnetic field gets stronger. 

In Figure H] the profile of the longitudinal (parallel to the magnetic field) 
and transverse (perpendicular to the field) components of the quark's electric 
current in a background of a specially-prepared smooth instanton-like {Q = 1) 
configuration subjected to an external magnetic field. The CME is character- 
ized by enhancement of the longitudinal current with respect to the transverse 
one. This property is clearly seen in the case of the instanton, Figure HI 

Coming back to real quantum configurations of the gluonic fields, we plot in 
Figure [5] (left) the expectation values of the fluctuations of each component of 
the electric current at zero temperature. We do not distinguish the global topo- 
logical charge of each conflguration since the CME effect appears also locally 
(an extended discussion of this point is given in Ref. [26]). In agreement with 
the CME features, we observe the dominance of the longitudinal components 
of the electric current with respect to the transverse ones. 

The CME crucially depends on the (local) imbalance between the left and 
right chiral modes given by the local chiral charge (x) = tp (x) 75^^ (x). The 
local strength of the chiral imbalance is characterized the expectation value 
of the chirality fluctuations, {pi), which is shown in Figure [5] (right) for three 
different temperatures. At T = the fluctuations quickly grow with the increase 
of the magnetic field. As the temperature increases the growth rate gets smaller, 
and in the deconfinement phase the rate is almost zero. In the confinement 
phase (T = 0) the chirality fiuctuations at strong enough magnetic fields are 
as large as the fluctuations in the deconflnement phase. This fact allows us to 
suggest that the CME may be observed in a cold nuclear matter as well. 
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qB =0.7 GgY^ ^5 = 1.8GeV^ 




Figure 3: A visual evidence of the CME at a typical T = gluon configuration: the 
excess of the positive (red) and negative (blue) electric charge due to the magnetic 
field (directed vertically) qB = 0.7 GeV^ (left) and qB = 1.8 GeV^ (right), Ref. 




Figure 4: An evidence of the CME at an instanton-like configuration of unit topolo- 
gical charge: a magnetic field induces the electric current j. The upper (green) 
and lower (blue) surfaces are the longitudinal (j| = + jf) and the transverse 
= il + j|) currents in the 12-plane (left) and in the 30-plane (right), Ref. [26] . 




Figure 5: (left) The squares of the transverse (ji, J2) and longitudinal (^3, jo) compo- 
nents of the current vs the magnetic field qB at T = 0. (right) The chirality squared 
vs the magnetic field qB at various temperatures T. Figures are from Ref. [26j . 
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5 Quark electric dipole moment 



The CME induces the electric dipole moment due to the spatial charge sep- 
aration along the direction of the magnetic field. Another manifestation of a 
nontrivial topological structure of the QCD vacuum is the appearance of the 
quark electric dipole moment directed along the axis of the external magnetic 
field. As in the case of the CME, the bulk average of the quark electric dipole 
moment is zero ((^'it^^') = 0) due to the global CP-invariance of the QCD 
vacuum. However, this anomalous effect may be seen on the event-by-event 
basis: locally the electric dipole moment of a quark may be large, while its sign 
may be alternating as the quark travels trough local fluctuations of the topo- 
logical charge. In order to evaluate the magnitude of the local electric dipole 
moment of the quark, we study the following connected expectation values: 

ii4f)iR = ((^3-(^^))')^^^-((^^-(^^))')^^^^„, i = E,M, (17) 

where crf{x) = tp (x) SjoV' (^) is the local density of the electric dipole moment 
[the density of its magnetic counterpart is (yf\x) = ^Eijki^ (x) T^jkip {x)]. 
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Figure 6: (left) The fluctuations of the longitudinal components the magnetic and 
electric dipole densities vs qB at T = and T = 1.12 T^. (right) The correlation of 
the electric dipole moment of quark (EDM) with the chiral density vs qB. 

In Figure [6] (left) we show values of the fluctuations of the longitudinal (i.e., 
directed along the magnetic field) component of the dipole moments (|17p . The 
fluctuations of the electric and magnetic dipole moments are equally strong in 
both phases. The transverse fluctuations of both electric and magnetic dipole 
moments are zero within our error bars. 

The electric dipole moment of the quark is closely related to the local chi- 
rality. In Figure [6] (right) we plot the (normalized) correlator of the electric 
dipole moment with the chiral density, c(/35,(j|') = {p5 cri) /[{pi) {{cri')'^)]^^'^ ■ 
At strong magnetic field we observe the very strong correlation of these two 
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quantities in both confinement and deconfinement phases, while at weaker mag- 
netic fields the correlation decreases. We also observed that the thermal fluctu- 
ations suppress the correlation function c [p5,a^), and that magnetic moment 
is not correlated with the local chirality. 

Thus, we found an evidence that the external magnetic field forces the quark 
to develop a local electric dipole moment along the direction of the field. 
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